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Abstract. We review all the cosmic shear results obtained so far, with 
a critical discussion of the present strengths and weaknesses. We discuss 
the future prospects and the role cosmic shear could play in a precision 
fS| . cosmology era. 
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r~| . 1. Introduction 

Or 

2 . The observation of gravitational lensing by large scale structures is a direct probe 

| of the matter distribution in the Universe. This method gives the most unbiased 

c3 picture of the matter distribution at low redshift compared to other techniques 

like cosmic velocity fields, galaxy distribution or Lyman-a forest studies. Indeed, 
these techniques rely on assumptions either like the dynamical stage of the struc- 
rS \ ture involved, or the properties of visible material versus dark matter biasing, or 

■ suffer of a poor sampling, or a combinaison of those. On the other hand, lensing 

by large scale structures suffers from practical difficulties, like its sensitivity to 
non-linear power spectrum predictions, or to the Point Spread Function correc- 
tions, which we will discuss later. In this review, we intend to give a present 
day picture of the cosmic shear research and to discuss the technical issues that 
could be a limitation. These technical limitations will certainly be overcome 
sooner or later, this is why a discussion of the role of cosmic shear for precision 
cosmology is also of interest. Although this paper is supposed to review the 
topic, there are already more than hundreds of publications on the cosmic shear 
subject alone. It is therefore difficult to address all aspects in details, and to 
mention everything (theory, simulations and observations). Instead, we choose 
to focus on observations, data analysis and related cosmological interpretations. 
By cosmic shear, we mean distorsion of the distant galaxies only. The magni- 
fication aspects of gravitational lensing by large scale structures, which is only 
at its beginning in terms of intensive observations, will not be reviewed. We 
apologize whose those of which work will not be discussed. 
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Figure 1. A light bundle and two of its rays C and CJ . £(w) is the 
physical diameter distance, which separates the two rays on the sky, 
viewed from the observer (w = 0). 

2. Linking galaxy shapes to theory 
2.1. Lensing by large scale structures 

Light propagation in the inhomogeneous universe We first have to define the 
homogeneous background universe notations (identical to Schneider et al. 1998). 
The metric of the homogeneous Universe is written in the form 



ds 2 = c 2 dt 2 - a 2 {t) dw 2 + fl{w)duj 



(1) 



where a(t) = (1 + z)^ 1 is the cosmic scale factor normalized to unity today, w(z) 
is the radial coordinate, and fx(w) is the comoving angular diameter distance 
out to a distance w(z). The radial distance w(z) is given by the redshift integral: 



f z c c f z 

w{z) = i dz 'n = h~ Jo y/(i + z')*n Q 



dz' 



+ (1 + z') 2 (1-^o-^a) + ^a 



■ (2) 



where Hq is today's Hubble constant, and the angular diameter distance /kt(u>) 
reads 



f K- 1 / 2 sin( v / ^u;) for K > , 

fx{w) = I w for K = , 

I {-K)- 1 ' 2 wah(y/=Kw) for K < , 

where K is the curvature 



(3) 



K 



{n + n A - 1) , 



(4) 
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Figure 2. Effect of n, 7 or uj on the displacement of two test particles 
1 and 2 located on a test ring (dot-dashed circle) with coordinates 
(d0,O) and (O,td0). 



with f^o and the mean density parameter and the vacuum energy today. 

Consider two light rays C and £' coming from a distant source and converg- 
ing to an observer, and define d9 as the observed angular vector between the two 
rays (Figur^). We use the Cartesian complex coordinates, so d6 = (d#i,id#2)- 
In the absence of any inhomogeneities along the line of sight, the physical dis- 
tance between the two rays at an angular distance fxiws) from the observer 
to the source is defined as £ = fxiws) d#. Due to the inhomogeneities, like 
clusters of galaxies, voids and filaments, the physical distance £ deviates from 
this simple relation, and can be linearized as: 

i = f K {ws)AdB = f K (w s )( K +J ~%) de 

= f K (w s )(K-iuj)de + f K (wshde\ (5) 

The matrix A is by definition the amplification matrix. The geometrical origin 
of this expression is easily understood when drawing how d#i and d#2 change 
with a small (but non- vanishing) k, 7 or uj (see FigureOJ). They are just numbers 
which describe the infinitesimal relative displacement of two rays C and £' . 

<*« « ( ) : ^ * ( -iX ) ; ^ K ( -tX ) < 6 > 

As we shall see now, they are the quantities which contain the cosmological 
information. By definition, k is called the convergence field, 7 the shear field, 
and uj the rotation field. When the shear is not expressed in the eigenspace 
(which is the case in Figure |2 for instance), 7 is a complex vector in general (see 
Figure 01). 

A light beam is a congruence of null geodesies, which are marked with re- 
spect to a fiducial (reference) geodesic having a tangent vector kn- The rays C 
and £ are two geodesies of the congruence, whose the separation £ = £1 + 1^2 
is defined as a space-like vector perpendicular to the wave- vector k^. As above, 
for an infinitesimal displacement along the congruence it is always possible to 
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decompose the geometrical deformation of the ray bundle into a uniform expan- 
sion 0, a shear a and a rotation W. This defines the well known optic scalars 
(Sachs 1961): 



where k^ u is the covariant derivative of the wave-vector and (//; v) and [//; v\ 
denote the symmetric and antisymmetric permutation of indices respectively. 
The evolution of the optic scalars along the congruence is completely determined 
by the optical scalar equations which depend on the gravitational field (Sachs 
1961): 



(La Z 
^ + a Q = T = C mvp k^t a P. (8) 

Here, the geodesic is parametrized with dA = dw/{\ + z) 2 . R^ u and C^ au p are 
the Ricci and the Weyl tensors respectively. t a is the complex null tetrad (or 
Sachs tetrad) such that t a k a = and t a t a = 1. Note that the first equation in 
(|5|) is nothing else but the Raychaudhuri equation for null geodesies. 

For an infinitesimal displacement along the congruence, the separation £ 
transforms according to Eq.©: 

|j = (9 -zW)€ + *•€*. (9) 

Differentiating Eq. @ and substituting Eq. (JHJ) leads to the evolution equation of 
£ along the congruence as a function of the gravitational fields 1Z and T: 



d 2 i _ (K-Re{T) Um{T) 



d\ 2 V Um(F) K + Re{F) 



£■ (io) 



The final step is to calculate 1Z and J- from the Ricci and the Weyl tensors for a 
Newtonian gravitational potential Straightforward but lengthy calculations 
give: 

n = -^—^- F=-^—(dl<s>-dl<s> + 2id 1 d 2 <s>) 1 (11) 

where a(w) is the scale factor of the unperturbed background metric, and w 
the radial distance. Using a perturbative expansion for the amplification matrix 

Aij = Afj + A\f + ... and for the gravitational potential $ = + ^ + 
Eq. ljlfljl can be solved iteratively. The homogeneous universe case corresponds 
to Aij = Aij = 5ij and $ = 0. It is then easy to obtain the general first order 
solution for the amplification matrix in the direction 9: 



Aussois's Cosmic Shear Lecture 



5 




Figure 3. Illustration of the first order effect of cosmic shear on a cir- 
cular background galaxy of radius Rq. The convergence is an isotropic 
distortion of the image of the galaxy, while the shear is an anisotropic 
distortion. 



Ajip) = s, j + 4\e) = 5, - 4 I" dw fK{W 7?{ K{wl) ^UK{ w ')e^ 

C z Jo JK\W) 

(12) 

where ws is the position of the source. Eq. ((T2l is the basic lensing equation 
used to calculate the distortion and the magnification of distant sources. This 
result is a first order expression and is only valid in the realm of the Born 
approximation where the lensing properties are calculated along the unperturbed 
light path (of direction 9). Therefore, all contribution coming from the lens- 
lens coupling are neglected. For most practical applications this is however an 
excellent approximation (Bernardeau et al. 1997, Schneider et al. 1998), as we 
shall see later. 

Back to the lensing effects (Eq|5J), the geometrical deformation of a light 
bundle can be expressed as an integrated effect along the line-of-sight: 

k=1 + \tt(A^); 7 = ^(4 1 i ) -4 1 2 ) + 2i4 1 2 ) ); w = 0. (13) 

These expressions show that a scalar perturbation will never induce a rotation 
of the light bundle at the first order (u = 0). Figure El shows the effect of cosmic 
shear on a distant circular galaxy, at the first order (rc <C 1 and 7 -C 1). 
It shows that the shear can be obtained from the measurement of the shape of 
galaxies. The practical methods to do this measurement will be discussed in 
Section 2.2. 
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Mean fields The second order derivatives of the gravitational potential field 
can be written as function of the mass density contrast 5, using the Poisson 
equation: 

V 2 d> = ^^5. (14) 
2a 

From Ea (|12|) . we get the convergence k{9) in the direction 9, as function of 5, 
integrated along the line of sight: 

3 H 2 n»s , f K (w-w')f K (w')5(f K (w')8,w>) 

k(9,w) = - — fl / dw — — — , (15) 

2 c Jo JK{w) a(w') 

with similar (but not identical) expressions for j(9). The sources have been 
assumed to be at a single 'redshift' w$, but similar expressions can be easily 
generalized for a more realistic redshift distribution. In that case, the lens- 
ing fields are integrated along the redshift with the proper source distribution 
p w (w)dw from to the horizon wn'- 

K ( e ) = o — n o / dwg(w)f K (w) , (16) 

2 V c / Jo avw) 



with 



9{w) 



I w 



i / , is fK{w - W) 

&w Pw(w fin - ( 17 ) 



Limber equation and small angle approximation We are primarily interested in 
the statistical properties of the lensing fields, which are given by the moments of 
the field. The variance is the first non trivial moment; its evolution with angular 
scale depends on cosmological parameters and on the geometrical properties 
of the Universe due to the light rays propagation. The mass density power 
spectrum Pzo{k) is defined as 

<5(k)5*(k')) = (27rf8 D (k-k')P 3D (k,w). (18) 
Likewise, one can define the convergence power spectrum P K (s): 

(«(s)«*(s')) = (2vr) 2 5 D (s - s') P K (s). (19) 

The time dependence in Ea (|18|) stands for the growth of structures. For an 
EdS Universe, it can be factorized, but in the general case it is more com- 
plicated, in particular in the non-linear regime where time dependence and 
scales are coupled. The jump from the 3-D wave vector k to the 2-D angu- 
lar wave vector s is ensured from the line of sight integration using the Lim- 
ber approximation (Limber, 1954). To simplify Ea (|16|) . it can be written as 
= / ^ q(w) 5(fK(w)9,w). In real space, the convergence correlation 
function £ K (A0) = (k(9)k(9 + A9)) can be eventually computed (Kaiser 1998): 

(k(0)k(0 + A0)) = J dw q(w) J dw' q(w') (5(f K (w)9,w)5(f K {w'){9 + A9),w')) 



Aussois's Cosmic Shear Lecture 



7 



~ J dwq 2 (w) J Aw' (5(f K (w)0,w)5(f K (w')(9 + A8), w')), (20) 

assuming that the selection function q(w) does not vary across the largest fluc- 
tuations of the density and that the fluctuations are much smaller than the 
distance of the sources. In order to express all cosmic shear 2-points statistics, 
we are in fact interested in the convergence power spectrum P K (s): 

P K (s) = J d6UO)e- is . (21) 
The density contrast 5(fK(w)0,w) = 6(r) can be expressed in Fourier space: 



= / (03 e ' lk± 9 SK{W) ^ D i +) ^ ^ (k) ' (22) 

where D^~\w) is the linear structure growth factor (see the next section non- 
linear power spectrum), and k = (kj_, £3), k^ is the wave-vector perpendicular 
to the line of sight. From this equation and Ea(|18j). one can express the density 
correlation function appearing in Eg 1)20(1 : 

(6(r)5*(r')) = J dk e^ 6 ^±< e + Ae ) /*(«') 

x e -i k 3 (w-w>) d[ + \w)D^(w') P 3D (k). (23) 

When, as in our case, the small angle approximation is valid ( |A0| < 1 — 
2 degrees), the transverse wave- vector kj_ carries most of the power at |k|; that 
is P3D(k) — P3D(k±) (Peebles 1980). The k 3 integration then gives a Dirac 
delta function 8d(w — w'). If we perform the variable change kj_ Jk{w) = s the 
convergence power spectrum becomes: 

PM-fto^M+H^Ftoi^. (24) 

J IkM 1 J \fK(w)J 

Back to the notations of Eq (|16j) . the convergence power spectrum finally writes 

P.WJffl'W-,,^/ • ;,). (25) 

4 V c J Jo a z (w) \Jk\w) J 

The shear power spectrum Py(s) is identical to this expression. The reason is 
that, in Fourier space, the quantities (k 2 ) and (I7I 2 ) are identical. This can 
be derived easily from Eq ((12() and Eq ((13() . with the derivatives replaced by 
powers in s's in Fourier space. As we shall see, this allows us to extract the 
convergence 2-points statistics directly from the data. Higher order statistics is 
a more difficult issue which will be discussed later. 
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Figure 4. The left panel is a 3-dimensional mass power spectrum 
for the linear (dashed) and non-linear (solid, using Smith et al. 2002) 
regimes when baryons are included. A value of = 0.05 was used. The 
right panel shows the induced convergence power spectrum (Eg 125(1 for 
the two dynamical regimes. Other parameters are O c dm = 0.25, = 
0.7,(78 = 0.9, h = 0.7, z sourcc = 0.8. 



Non-linear power spectrum The normalization of the mass density power spec- 
trum P^£] is defined in the conventional way, by computing the mass density 
variance within a sphere of 8 h Mpc radius at redshift zero: 

4 = (4) = ^3 / dkP 3D (k,0)\W(kR)\ 2 , (26) 

where W(kR) = ( ^H^TT^ ~~ cos(kR)^j is the Fourier transform of the top- 

hat window function of radius R. The transition from the linear to the non-linear 
scales is identified by <7g ~ 1. In the linear regime, where the density contrast 
of the mass distribution is low (5 <C 1), the fluid equations describing the 
structure growth can be solved perturbatively, and one obtains for the growing 
mode: 



P 3D (k,w) = D[ + \w) ' P 3D (k) 



(27) 



with, 



^ Jo 



da 



a 3 H(a) ' 



(28) 



In the non-linear regime, the structure growth cannot be solved analytically and 
its description must rely on non-linear models (Peacock & Dodds, 1996, Smith et 
al. 2002), following an original idea of Hamilton et al. (1991). Non-linear predic- 
tions of the matter power spectrum are performed from the knowledge of the spa- 
tial 2-points correlation function of the galaxies £2(1") = tt^3 I dk P(k) e~* k r . 



An accurate measurement of ^2( r ) is given for instance by the 2dF (Percival et 
al. 2001) or the SDSS surveys (Dodelson, S., et al. 2002): 
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^)=(^) 7 , (29) 

with ro = 4.3±0.3 h M.pc and 7 = 1.71±0.06. The stable clustering hypothesis 
stipulates that at very small scales (strong non- linear regime) , the internal profile 
of clusters of galaxies remain constant with time for any cosmological model, and 
that the cluster distribution is driven by the cosmic expansion. This means that 
the correlation function is fixed in proper coordinates, but its amplitude evolves 
as a volume effect like (1 + z) -3 . At large scale (linear regime), the correlation 
function follows the perturbation theory. Since the correlation function ^( r ) 
behaves like r~ 7 for any cosmological model, we therefore have the two following 
limiting cases (Peacock 1999): 

^2(^5 z) oc (1 + z) J (l + z)~ 3 non — linear (30) 



r (+) 1 2 

^{r, z) cc D\ (w) linear (31) 

A mapping from the linear to the non-linear scale has been conjectured (Hamil- 
ton et al. 1991, Peacock & Dodds 1996, Smith et al. 2002), and calibrated using 
N-body simulation. The transition from linear to non-linear scales is described 
by a few slowly varying functions that depend on cosmological parameters. The 
same argument applies to the 3-D power spectrum, which is needed for cosmic 
shear predictions down to small scales fEa 125(1 (Peacock & Dodds 1994). Figure 
0] is an example of 3-dimensional and convergence power spectra in the linear 
and non-linear regimes. A fair amount of baryons was included (using CAMB, 
Lewis et al. 2002), in order to show that the baryon oscillations, which are clearly 
visible on the 3D spectrum, are severely diluted in the projected spectrum. 

2-points statistics In practice, the variance of the convergence (or shear, which 
is the same) is computed within a given smoothing window U{9) of radius C , 
which can be written: 



(* 2 )e c = <(/ dV U{6>)k{6>) 

d 2 0' U(0') I d 2 tf U{&)(k{&)k{*)). (32) 



If we express the convergence from its Fourier transform k(8) = J d 2 s k(s) e l ^' s 
and using Ea (|19|) . we obtain: 



(A 



/ dV U(9>) J d 2 tf UQ) J ^ e^ 9 '-^P K (s) 



2tt 



ds sPJs) 



d$$U($)J (s$) 



(33) 
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Figure 5. In order to compute the shear variances, the galaxy ellip- 
ticities are smoothed within a window (dashed red) of fixed radius 9 C 
(left). The shear variance will show up as an excess of galaxy align- 
ment with respect to random orientation. The right panel shows the 
profile of the two filters one usually consider, top-hat (solid line) and 
compensated (dashed line). On the left, the axis (et,e r ) correspond to 
the local frame attached to each individual galaxy, on which the galaxy 
ellipticity components can be projected out to give an estimate of the 
tangential 7$ and radial shear 7 r . 



This expression is general, and can be applied to any smoothing window U(9). 
Since P-y(s) = P K (s), it also expresses the shear variance ( , y 2 )e c - As illustrated 
in Figure [5] , we are primarily interested in a top-hat filtering, for which, 



(7 2 



/ ds s P K (s) 



Ji(s9 c ) 

S0 a 



(34) 



and in the compensated filtering having J C &9 6 U{9) = (zero mean). The 
choice of U{9) is arbitrary, provided it has a zero mean. Here we use the expres- 
sion (Schneider et al. 1998): 



(35) 




so the variance of the convergence with this filter is: 



{Ml p ) = —J dssP K (s) 



J 4 (s9 c ) 



(36) 



The nice feature of the compensated filter is that it is a pass-band filter, which 
means that the variance Eq ()36|) is a direct estimate of the convergence power 
spectrum in real space. Note that the power is estimated around s ~ 5/# c . 
Furthermore, it can be estimated directly from the ellipticity of the galaxies, 
without a reconstruction of the convergence field. This remarkable property has 
been demonstrated by Kaiser et al. (1994), who have shown that Ea (|36|) can be 
obtained from a smoothing of the tangential component of the shear field 74 : 
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0.1 1.0 10.0 100.0 

sO 

c 

Figure 6. Top-hat (dashed line) and compensated (solid line) filters 
in Fourier space. This plot illustrates the fact that the compensated 
filter is a pass-band filter, and therefore is a broad-band estimates of 
the convergence power spectrum in real space. 



M ap = / tie Q(8) 7t , (37) 
Jo 

where 

Q(d) = ^lJ o 9c ae'9'u(9')-u(e). (38) 

The tangential shear 74 can be obtained from the projection of the galaxy ellip- 
ticity on the local frame (Figure EJ). 

Another 2-points statistics of interest is the shear correlation function (7 • 
7)e c . It consists in calculating the sum of the shear product of all possible 
galaxy pairs separated by a distance C . Using the shear field version (i.e. for 
7) of Ea (|15|) . one can show that (Blandford et al. 1991, Miralda-Escude 1991, 
Kaiser 1992): 

h ■ f)e c = ^J dss P -( s ) MsOc). (39) 

One can also compute the shear correlation functions of the projected compo- 
nents of the shear, (74 74), (7,. 7 r ). For symmetry reasons (74 7,.) = 0. On 
the other hand, the two former correlation functions are not equal, because the 
gravitational shear is generated by a scalar potential, implying that the projec- 
tions on the local frame of the shear components are not equivalent. We can 
show that: 
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Figure 7. Lensing statistics predictions for the cosmological model 
used in Figure |1] Both linear (dashed) and non-linear (solid lines) 
regimes are represented. On the bottom-right plot, the thick dashed 
and solid lines are the full shear correlation function. The cosmological 
model is the same as Figure HJ 



(it lt)e c = Jds s P K (s) [J (s9 c ) + J 4 (s6 c )] 

(7r7rk = Jds s P K (s) [J (se c ) - J 4 (s0 c )] (40) 

One usually denotes £+(0 c ) = (it It) + (lr lr), and £_(0 C ) = (7* 7*) - {jr lr)- 
We have, of course, £+(0 c ) = (7 • l)e c - 

Figure [7| shows the linear and non-linear predictions for all the statistics 
defined here, for a particular cosmological model. 

Dependence on cosmological parameters It is obvious from Eq ()15|) . Ea (|25|) and 
Ea(|26|l that the cosmic shear signal depends primarily on the source redshift ws, 
then on the mean density parameter Qq, and on the slope and the normalization 
(as) of the mass power spectrum. To explore the parameter dependence of the 
cosmic shear signal, we assume the Cold Dark Matter model, with a power spec- 
trum parameterized with the slope parameter T. We allow the four parameters 
(Qo, z s , r, as) to vary, and we compute the likelihood C(Qq, z s , T, as \ d) of the 
parameters knowing the data d. The data vector is for instance the aperture 
mass or any other statistic: 
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Figure 8. 1-cr, 2-cr and 3-cr confidence contours for the maximum 
likelihood analysis on the four parameters f2 m , erg, T and the source 
redshift parameter z s (see text). The six possible pairs of parame- 
ters are displayed. On each figure, the two hidden parameters are 
marginalized such that Q m £ [0.2,0.4], cr 8 £ [0.8, 1.1], T £ [0.1,0.3] and 
z s £ [0.4,0.5], and the cosmological constant is fixed to = 1 — f2 m . 
The reference model is O m = 0.3, as = 1, T = 0.21 and z s = 0.44. The 
survey area is A = 16deg 2 , the galaxy ellipticity r.m.s. is 0.3, and the 
correlation functions are measured in the range 0'6 < 9 < 30'. 



C 



(2vr)«/ 2 |S| 1/2 



exp 



(d 



s) T S- 1 (d-s) 



(41) 



where s is the fiducial model vector and S := ((d — s) T (d — s)) is the covari- 
ance matrix. Figure |H1 and show the parameter dependence one expects for a 
survey covering 16 square degrees up to the limiting magnitude Iab = 24, for 
two different choices of priors. The signal also depends on other cosmological 
parameters (f^, Q\, ft u ,...), albeit to a lower extend. For precision cosmology, 
all parameters are relevant, but the first constraints obtained so far from cosmic 
shear are on the main four parameters (Qq, z s , T, eg). 
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Figure 9. Same as figure |Hl with strong priors: in each figure, the two 
hidden parameters as assumed to be known perfectly. These plots show 
the degeneracy directions among all the possible pairs of parameters 
obtained from S7 m , erg, T and z s . 



2.2. Galaxy ellipticities and estimators 

Ellipticity of the galaxies As mentioned in the previous Section, the cosmic 
shear signal is measured from the shape of the distant lensed galaxies. It is 
quantified from the ellipticity e. The raw ellipticity e of a galaxy is measured 
from the second moments 1™ of the surface brightness /(#): 

6 = (%K^f 5 ^)) ' hj = I 

The window function W{9) suppresses the noise at large distances from the 
object center. The cosmic shear signal can also be measured using gravitational 
magnification from the relative size and number count of the lensed galaxies, but 
this is out of the scope of this paper. Here, we only focus on the gravitational 
distortion effect. If one could measure the shape of the galaxies (with W{6) = 1) 
perfectly without any systematics coming from the telescope tracking and the 
optical defects, and if the galaxies were only lensed, then the observed ellipticity 
would be related to the source ellipticity as 
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Figure 10. Right plots: simulated cores of anisotropic PSF's. Left 
plots: simulated outer part of the PSF (Erben et al. 2001). The 
different PSFs are computed from ray-tracing through the telescope 
optic. 



pSt/urue _i_ p. 

e obs = _^ (43) 

]_ _|_ ^source . g 

where g = 7/(1 — k) is the reduced shear, and e obs is the observed ellipticity, 
^source j g source (unobserved) ellipticity. For nearly all cosmic shear applica- 
tion, the lens fields are small (|g|,K <C 1) and the linear approximation is valid 

e obs ^ e source _j_ 

Unfortunately, the ellipticity of the galaxies measured on the images are 
contaminated by atmospheric and instrumental distortions of the Point Spread 
Function (PSF) that also produce coherent non-gravitational elongation pat- 
terns, even on stars. Such example of PSF is displayed on Figure 1101 and the 
measured coherence of the PSF distortion on a real field is shown on Figure 
1111 This is a critical issue, for instance there were two early tentatives to mea- 
sure the gravitational lensing by large scale structures, which failed because the 
image quality was very low and the PSF correction not accurate (Valdes et al. 
1983, Mould et al. 1994). Since then, various methods have been developped to 
correct for the non gravitational source of galaxy alignment, which followed the 
improvement of image quality: 

• Kaiser et al. (1995), a method which treats the PSF convolution analyti- 
cally to the first order. It is called KSB. 

• Bonnet & Mellier (1995), which combines galaxy image simulation and 
optimal weighting of the isophotes. 
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• The auto-correlation function (Van Waerbeke et al. 1997), similar to Bon- 
net & Mellier (1995), but applied to the auto-correlation of the image of 
the galaxies to avoid some problems associated with the galaxies. 

• Kuijken (1999), a method which parametrizes the PSF and the galaxies 
with analytical functions, and try to match the convoled profile to the 
data. 

• Kaiser (2000), extended KSB, which circularises the PSF before the isotropic 
correction. 

• Modified KSB (Rhodes et al. 2000), is the KSB approach, applied on the 
galaxy moments instead of the ellipticities. 

• Bernstein & Jarvis (2002), is first a circularisation of the PSF, and then the 
convolved profile is analysed using a reduced set of orthogonal functions 
(Laguerre polynomials). 

• The shapelets approach (Chang & Refregier, 2002), is a kind of Principal 
Components Analysis, using orthogonal Hermite polynomials functions to 
decompose the convolved galaxy images (see also Bertin 2001 for a PC A 
approach) . 

The most popular, and certainly the most intensively tested 1 (Erben et 
al. 2001, Bacon et al. 2001), is the KSB approach. It is a very simple and 
powerful correction based on the first order effect of a convolution. The idea is 
that we can write the first order effect of the shear and of the PSF convolution 
analytically as: 

e obs = source + . ^ + psm . ^ (44) 

where P 7 and P sm are tensors computed on the image (see Kaiser et al., 1995), 
e* is the star ellipticity at the galaxy location, and 7 is the shear signal we want 
to measure. Assuming that the galaxies are isotropically oriented in the source 
plane, we have (e source ) = (which is valid even if the galaxies are intrinsically 
correlated), therefore the shear estimate from the measured galaxy ellipticity is 
given by: 

7 = p- 1 (e o6s - P sm • e*) . (45) 

We discussed in the previous section how the shear (7) could be splitted into 
a radial and a tangential component 7 r and 74 when projected onto the local 
frame of the aperture (Figure |SJ). Figure ^] shows the relation between the 
components e = {e\,e<i) of a galaxy, and its orientation. If we identify (ei,e2) 
to (et, e r ), we obtain the orientation in the local frame. 



*A realistic image simulation software is available at http:/ /affix. iap.fr/soft/skymaker/index. html 
and a realistic catalogue generation at http://affix.iap.fr/soft/stuff/index.html 
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Figure 11. Uncorrected (left) 
one of cosmic shear fields. 
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and corrected(right) star ellipticities in 



E and B modes The gravitational field is supposed to be completely dominated 
by a scalar gravitational potential at low redshift. The consequence is that 
only curl free modes for the shear are allowed. Any significant curl component 
should be interpreted as a (bad) sign of residual systematics in the data. Figure 
EI] shows the E mode generated by over-densities (top-left) and under-densities 
(top-right). The two bottom curl modes are not allowed. Using the statistical 
properties of these patterns and the (et,e r ) conversion from Figure IT21 it can 
be shown that the E modes correspond to the aperture mass (Mf p ), and the B 
mode to the aperture mass with the galaxies 45 degrees rotated (such rotation 
corresponds to a switch e t — ► e r ; e r — ► — ej). This is easy to understand: if there 
is no B mode, then switching the E into B, and B into E modes kills the signal 
measured with the aperture mass statistics. 

Aperture mass from the shear correlation function Because the E/B mode sep- 
aration provides a direct and robust check of systematics error residuals, it is 
widely believed to be the most reliable statistics. In order to compute it, there is 
fortunately no need to draw a compensated filter across the data and to average 
the shear variance; otherwise, this could be terribly complicated with real data 
because of the complex shape of the masks (see Figure HI)) . Variances and cor- 
relation functions can be expressed one into another (since they are only linear 
combinaison one to another). The E mode aperture mass is given by 

/•20c [2e c 
{M% p )=irj rdrW{r)5, + (r) + vr J rdr>V(r)f_(r), (46) 

where W(r) and VV(r) are given in Crittenden et al. 2002 and Pen et al. 2002. 
The i?-mode is obtained by changing the sign of the second term in Ea. (|46|) . The 
correlation functions and £_(r) are compute from the tangential and radial 

correlation functions (see Eal40|). In order to estimate the shear correlation 
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e r 




Figure 12. Value of (ej,e r ), or (e±, e-i) in Cartesian coordinates, as 
a function of the shape of a galaxy with respect to the local frame at- 
tached to the galaxy. Note that the ellipticity is invariant by a rotation 
of 7r, and not 2tt, this is why ej < and e r = for a vertical galaxy for 
instance. 



functions, let 8i be location of the i-th galaxy, its ellipticity e(0j) = (ei,e2), 
and the weight Wi. The ellipticity is an unbiased estimate of the shear ~f(0i). 
The quantity measured from the data are the binned tangential and radial shear 
correlation functions. They are given by a sum over galaxy pairs (8i,0j) 



WiWje t {Oi) ■ e t (9j) Y w iWje r {6i) ■ e r (6j) 
&(r) = ^ = ; Mr) = ^ = , (47) 



where r = \0i — 0j\, and (ej, e r ) are the tangential and radial ellipticities defined 
in the frame of the line connecting a pair of galaxies. The weights wi are usually 
computed for each galaxy from the intrinsic ellipticity variance o\ and the r.m.s. 
of the ellipticity PSF correction of. For example, van Waerbeke et al (2000) 
measured a e ~ 0.4 from their CFHT data, and defined the weights as: 



Wi 



at + a, 



2 " 



(48) 



To compute a t for each galaxy, the galaxy size-magnitude parameter space is 
divided into cells of constant object number (typically 30 galaxies per cell). For 
each cell the r.m.s. of the ellipticity correction among the galaxies in the cell 
is computed. This choice of parameter space is motivated by the fact that the 
isotropic PSF correction (the P 7 term in Eg 1441) is mainly sensitive to the size 
and magnitude of the galaxies. 
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Figure 13. Top patterns: shear curl free modes (E modes) allowed 
by gravitational lensing. Bottom patterns: curl modes (B modes) not 
allowed from a scalar gravitational potential. Only the E modes gives 
the signal of the aperture mass statistics {M^ ). 

3. 2-pts statistics 
3.1. Measurements 

There are now several evidences of the cosmological origin of the measured signal: 

(a) The consistency of the shear excess variance measured from different 
telescopes, at different depths and with different filters. This is summarized 
on Figure Hoi The first detections were obtained by Bacon et al. 2000, Kaiser 
et al. 2000, Van Waerbeke et al. 2000, Wittman et al. 2000. Since then, 
several measurements have been done in different observing conditions, which 
are summarized in Tabled 

(b) On a single survey, the self consistency of the different types of lensing 
statistics given by Eas. (|:UMK9l4ol) . This was done on the VIRMOS-DESCART 
survey 2 , and it is shown in Figure ITol (Van Waerbeke et al. 2001). 

(c) The comparison of the E and B modes measurements (to higher accu- 
racy than in (b)) between a deep and shallow survey for the VIRMOS-DESCART 
and RCS 3 surveys (Van Waerbeke et al. 2002, Hoekstra et al. 2002). This is 
shown on Figure IT7I More recently, the E and B modes have been also measured 
in other surveys (Brown et al. 2003, Jarvis et al. 2003, Hamana et al. 2003), 
which supports the cosmological origin of the signal, showing also the already 



2 http://www. astrsp-mrs.fr and http://terapix.iap.fr/DESCART 
3 http: //www. astro. utoronto.ca/ gladders/RCS/ 
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Figure 14. Mask area for real data. These holes of various sizes make 
the mass reconstruction very challenging. Each CCD chip is about 
T x 14'. Entire CCD's had to be removed because of bright stars and 
residual fringes patterns. 



small amount of residual systematics achieved with today's technology. The E 
and B mode measurements should now be considered as the most robust proof 
of the cosmological origin of the signal, and a quantitative test of systematics. 

(d) The lensing signal is expected to decrease for low redshift sources, as 
consequence of the lower efficiency of the the gravitational distortion. It cor- 
responds to a change in w s in EaQ15|). or equivalently a change in the mean 
source redshift with Ea(|25|). This decrease of the signal has been observed for 
the first time with the comparison of the E mode amplitude of the VIRMOS 
survey aperture mass (see Figure IT7j). which has a source mean redshift around 
0.9, to the RCS which has a source mean redshift around 0.6. The expected 
decrease in signal amplitude is about 2, which is what is observed. This is a 
direct evidence of the effect of changing the redshift of the sources, a kind of 
3-D cosmic shear effect. 

(e) Space images provide in principle a systematics- low environment, and 
even if the observed areas are still smaller than ground based observations, space 
data provide ideal calibrations of the cosmic shear signal (Rhodes et al. 2001, 
Haemmerle et al. 2002, Refregier et al. 2002), which are in agreement with 
ground based measurements (see Figure Hoi the HST points). 

3.2. Constraints 

The standard approach is to compute the likelihood of a set of n parameters 
(j>i,P2, ■■■,Pn), knowing the data vector d, as written in Ea lflTj) . As the data 
vector, it is natural to choose the aperture mass variance as a function of scale 
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Figure 15. Compilation of recent results of top-hat shear variance 
measurements from several groups (Courtesy Refregier et al. 2002). 



(M^ p ), because the signal is splitted into gravitational lensing and systematics 
channels (the E and B modes). The B mode measures an estimate of the 
contamination of the E mode by systematics. The E and B modes do not 
equally contribute to systematic, but we know, from the measurement of the 
modes on the stars, that they are very similar. If the B mode is not consistent 
with zero (which is the case for all surveys at the moment), it is important to deal 
with it properly when estimating the cosmological parameters. Unfortunately 
it is not yet clear what the best approach is: some groups (Van Waerbeke et 
al. 2002, Hoekstra et al. 2002, Hamana et al. 2003) added the B mode in 
quadrature to the E errors, taking into account the correlation between various 
scales. The B mode has been subtracted first from the E mode in Hoekstra et 
al. (2002), but not in Van Waerbeke et al. (2002). This might probably result in 
a slight bias for high <rg values in the later. Unfortunately we have no guarantee 
that the B subtraction is the right correction method. Recently Jarvis et al. 
(2003) marginalised the probabilities over E — B to E + B taken as the signal, 
which is more likely to include the 'true' B mode correction one has to apply. 

Figure UHl shows 4 the E and B modes that have been measured so far, using 
the aperture mass only (this is the only statistic which provides an unambiguous 
E and B separation, Pen et al. 2002). The two deepest surveys have large scale 
B mode contamination (Van Waerbeke et al. 2002, Hamana et al. 2003), and 
the two shallow surveys have small scale contamination (Hoekstra et al. 2002, 
Jarvis et al. 2003). 



4 The B mode peak at 10' in Hamana et al. (2003) is due to a PSF correction error over the 
mosaic. It is gone when the proper correction is applied, Hamana, private communication. 
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Figure ITUl shows the joint fi m , erg constraints obtained from the measure- 
ments of Figure 1171 They are obtained only when comparing the measured 
lensing signal to the non- linear predictions. Unfortunately, the actual surveys 
are not yet big enough to probe the linear scales accurately. The non-linear 
power can be computed numerically (Smith et al. 2002), but its precision is 
still uncertain. Recent investigations show that a 10% r.m.s. uncertainty is ex- 
pected, which means that the cosmological parameters cannot be known with 
better precision for the moment. According to the Figure [7J the transition scale 
between the linear and non-linear regimes is around 1 degree. The consequence 
is that the quoted mass normalization eg is sensitive to the validity of the non- 
linear mapping at small scale. In this respect, Jarvis et al. (2003) are less 
contaminated by this problem because they used the lensing signal from 30' to 
100' to constrain the mass normalization. 

Table ^ summarizes the ag measurements for all the lensing surveys pub- 
lished so far. For simplicity it is given for £l m = 0.3. Despite the differ- 
ences among the surveys, it is worth to note that the results are all consistent 
within 2.5a between the most extreme cases, when poorly known parameters 
are marginalised. 
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Figure 16. Measurement of all the 2-points statistics in the same 
survey, VIRMOS (Van Waerbeke et al. 2001). Top left: top-hat vari- 
ance. Top right: aperture mass E and B modes. Bottom left: full 
shear correlation function. Bottom right: projected shear correlation 
functions, survey. Right: E (top) and B (bottom) modes measured in 
the RCS survey. The B mode is low and the E mode compatible with 
the predictions for the aperture mass statistics. The lines are fiducial 
models which indicate the relative deviations between the statistics to 
the models. 
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Figure 17. Left: E (top) and i? (bottom) modes measured with all 
the most recently reduced data in the VIRMOS survey (Van Waerbeke 
et al. 2002). Right: E (top) and B (bottom) modes measured in the 
RCS survey (Hoekstra et al. 2002). The B mode is low and the E 
mode compatible with the predictions for the aperture mass statistics. 
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Figure 18. Plot showing the relative amplitude of the aperture mass 
E and B modes (points without and with error bars respectively) for 
all the surveys where the aperture mass has been measured (Hoekstra 
et al. 2002, Van Waerbeke et al. 2002, Jarvis et al. 2003 and Hamana 
et al. 2003) (picture taken from Jarvis et al. 2003 and extended). The 
result of Hoekstra et al. (2002) is for the full magnitude range, while 
in Figure IT71 right panel, it is for the galaxies used for the cosmic shear 
analysis. 
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Figure 19. The solid lines on each plot show the 1, 2 and 3<r contours 
of the VIRMOS and RCS survey, from the measurements shown in Fig- 
ure El The contours have been marginalised over the source redshift 
and the slope of the matter power spectrum as described elsewhere 
(Van Waerbeke et al. 2002, Hoekstra et al. 2002). 



Table 1. Constraints on the power spectrum normalization "as" for Q, m = 0.3 for a flat Universe, obtained from 
a given "statistic". "CosVar" tells us whether or not the cosmic variance has been included, "E/B" tells us whether 
or not a mode decomposition has been used in the likelihood analysis. Note that Van Waerbeke et al. (2001) and 
Brown et al. (2003) measured a small B-mode, which they didn't use in the parameter estimation. z s and T are the 
priors used for the different surveys identified with "ID". Note also the cosmic shear results obtained by Kaiser et 
al. (2000) and Haemmerle et al. (2002), which are not in the table here because they reported a shear detection, not 
a o"8 measurement. 



ID 




^8 


Statistic 


Field 


miim 


CosVar 


E/B 


Zs 


r 


Maoli et al. 01 




1.03 ±0.05 




VLT+CTIO 




no 


no 




0.21 








+WHT+CFHT 












Van Waerbeke et al. 


01 


0.88 ±0.11 


(r 4 ), m <M£> 


CFHT 8 sq.deg. 


1=24 


no 


no (yes) 


1.1 


0.21 


Rhodes et al. 01 




91+ lJ - 2b 


£(r) 


HST 0.05 sq.deg. 


1=26 


yes 


no 


0.9-1.1 


0.25 


Hockstra et al. 02 




0.81 ±0.08 


<7 2 > 


CFHT+CTIO 


R=24 


yes 


no 


0.55 


0.21 








24 sq.deg. 












Bacon et al. 03 




0.97 ±0.13 


£(r) 


Keck+WHT 


R=25 


yes 


no 


0.7-0.9 


0.21 








1.6 sq.deg. 












Refrcgier et al. 02 




0.94 ±0.17 


<7 a > 


HST 0.36 sq.deg. 


1=23.5 


yes 


no 


0.8-1.0 


0.21 


Van Waerbeke et al. 


02 


0.94 ±0.12 




CFHT 


1=24 


yes 


yes 


0.78-1.08 


0.1-0.4 








12 sq.deg. 












Hoekstra et al. 02 




o.9i±8:E 


<7 2 UM 


CFHT+CTIO 


R=24 


yes 


yes 


0.54-0.66 


0.05-0.5 








53 sq.deg. 












Brown ct al. 03 




0.74 ±0.09 


<7 2 >, 


ESO 1.25 sq.deg. 


R=25.5 


yes 


no (yes) 


0.8-0.9 




Hamana et al. 03 




(2a)0.69™ 


<M 2 p U(r) 


Subaru 2.1 sq.deg. 


R=26 


yes 


yes 


0.8-1.4 


0.1-0.4 


Jarvis et al. 03 




(2<7)0.71t££ 


< 7 2 ), £(r), <M 2 p > 


CTIO 75 sq.deg. 


R=23 


yes 


yes 


0.66 


0.15-0.5 
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4. 3-pts statistics 

So far, we only discussed the 2-points statistics, but recently higher order statis- 
tics have been also developed for cosmic shear (Bernardeau et al. 1997, Jain 
& Seljak 1997). If we were able to reconstruct the convergence from the shear 
(ellipticity) measured on the galaxies, one could for instance measure the top- 
hat smoothed higher order statistic easily. For instance, the skewness of the 
convergence, which is defined as 



is of great interest because this suited ratio of moments makes this statis- 
tic nearly independent of the normalization and shape of the power spectrum 
(Bernardeau et al. 1997). A pedagogical way to compare the second and third 
moments is to compute (k 2 ) and S3 (k) in the perturbation theory, and with a 
power law power spectrum. In that case, one finds 



These are only approximated relations, which are not valid in the real (non- 
linear) world, but it shows that the skewness provides a direct geometrical test 
(dependence on Qq), as long as we know the redshift of the sources z s . Combined 
with the second order moment, the degeneracy between the power spectrum 
normalization and the density parameter can be broken with the cosmic shear 
alone. 

The skewness can also be predicted in the non-linear regime, as for the 
2-points statistics, using a non-linear extension of the bispectrum (Scoccimarro 
et al. 2002, Van Waerbeke et al. 2002). The problem with the skewness of the 
convergence is that it cannot be measured on the data directly, and one needs to 
reconstruct k from the shear first. This process is unfortunately sensitive to the 
survey geometry because the projected mass reconstruction is essentially a non- 
linear process. Given the typical observed field geometry, as shown on Figure 
1141 it is yet impossible to perform a mass reconstruction with the accuracy 
required to measure the cosmic shear effect. One possibility for avoiding the 
mass reconstruction (that is try to make the map making a local process) is 
to compute the third moment of the aperture mass (Schneider et al. 1998). 
Unfortunately, in that well, the complicated survey geometry make it 

difficult to measure an accurate third moment (Mf ). 

The alternative is to measure a third moment on the shear field itself, but 
this cannot be done in a trivial way, since for evident symmetry reasons, any 
odd moment of the components of a vector field vanishes. One has to built 
explicitly non-trivial measures of the third moment of the shear, which has been 
recently proposed. So far, two of the proposed estimators lead to a measurement 
(Bernardeau et al. 2002 and Pen et al. 2002). 




(49) 




(50) 



(51) 
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Figure 20. Average shear pattern obtained around a galaxy pair lo- 
cated at (xi, X2). A third galaxy is located at x, its shear vector is 
projected along the vertical axis, it is called 7t(x). The shear 3-points 
function (7(xi) • 7^2)74 (x)) is averaged inside the ellipse indicated by 
the solid line. 



In Bernardeau et al. (2002) the idea is to identify regular shear patterns 
around any pair of lensed galaxies. A pair is identified by the two galaxy posi- 
tions xi and X2, and any location around the pair by x. For a fixed pair (xi, 
X2), we are interested in the average shear at x. 

Figure |2U] shows the typical shear pattern observed around a galaxy pair 
located at (xi, X2). Ray tracing simulations demonstrate the stability of this 
shear pattern, which is almost independent on the cosmological model and the 
pair separation. A natural 3-points function to calculate is the average of the 
product of the shear correlation function 7(xi) • 7(x2) with a projection of the 
shear of the third galaxy 7(x). It is obvious from Figure l2*Ul that the projection 
is optimal when performed along the vertical axis. For a fixed pair location (xi, 
X2), the 3-points function ^(x) is defined as: 

&(x) = <7(xi)-7(x 2 )7i(x)>, (52) 
and the quantity we measure is: 

6(|x 1 -x 2 |)=/ ^£ 3 (x'). (53) 
Jew. Veil 

Figure Ell shows the result on the VIRMOS-DESCART survey. The treat- 
ment of the B mode is still uncertain, and the redshift uncertainty still too large, 
which makes very difficult the interpretation in terms of cosmological parame- 
ters. However Figure shows that the order of magnitude, and the slope of 
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Figure 21. On the left, results for £ 3 (di 2 )/6(di2) 2 for the VIRMOS- 
DESCART survey (dot-dashed lines: E — B mode for the 2-points 
function, solid line: E + B mode for the 2-points function). This 
is compared to rCDM and OCDM results (dotted and dashed lines 
respectively). Right plot: dashed line is £3(^12) for the VIRMOS- 
DESCART survey, compared to the same quantity measured on the 
stars. 



the signal are consistent with the expectations. For instance, the signal from 
the stars before PSF correction is completely different in shape and amplitude. 

In Pen et al. (2002), the idea is to compute the convergence aperture mass 
3-points function from an integral of the shear 3-points function. This solution 
avoids the problem of drawing cells across a complex field geometry and presents 
the advantage to estimate the third moment of the convergence k, which is the 
field of physical interest. Unfortunately, its measurement is still very noisy, 
because it uses a compensated filter that removes the low frequency modes for 
any target frequency (which is not the case for a top-hat filtering) . The resulting 
skewness is shown on Figure l2*2l and is consistent with £Iq < 0.4 at the 90% level. 

Other approaches have been proposed (Zaldarriaga & Scoccimarro 2002, 
Schneider & Lombardi 2003, Takada & Jain 2003) which all deal with trying 
to optimize the signal-to-noise by looking for the best galaxies triangle config- 
urations containing the highest signal. They have not yet been applied to the 
data. 



5. Galaxy biasing 

A direct byproduct of cosmic shear observations is the measure of the light/mass 
relation, the so-called biasing parameter b defined as the ratio of the galaxy 
density contrast to the dark matter density contrast 

Sgai = b S mass . (54) 

This is in fact a highly simplified model, which assumes that the biasing does 
not vary with scale and redshift, and that the relation between mass and light 
is deterministic. While in the real world, none of these assumptions are correct, 
this model has the advantage to be tractable analytically, and to provide an 
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Figure 22. Skewness of the convergence as measured in Pen et al. 
(2002), on the VIRMOS-DESCART survey. The overall significance 
of the measurement if 3.3 a, which was computed using Monte Carlo 
sampling of the errors from ray-tracing simulations. A comparison of 
the signal with tese simulations shows that Qq < 0.4 at the 90% level. 
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Figure 23. Left plot: The measurements of (AT 2 ) (panel a), and 
(AAM ap ) (panel b) as a function of angular scale from the RCS data. 
Panel c shows (M ap ) as a function of angular scale from the VIRMOS- 
DESCART data. The error bars for (M ap ) have been increased to 
account for the unknown correction for the observed "B"-mode. For 
reference, a few models have been plotted, assuming 6=1 and r = 1, 
for an OCDM cosmology (dotted line; Q m = 0.3, Q\ = 0, as = 0.9, and 
T = 0.21) and a ACDM cosmology (dashed line; Vt m = 0.3, ft A = 0.7, 
o"8 = 0.9, and V = 0.21). Note that the points at different scales 
are only slightly correlated. Right plot: The measured value of the 
galaxy-mass cross correlation coefficient r as a function of scale for the 
ACDM cosmology, (b) The bias parameter b as a function of scale. 
The upper axis indicates the effective physical scale probed by the 
compensated filter at the median redshift of the lenses (z = 0.35). 



average biasing estimates, which is still very useful. Nevertheless, it is possible 
to go beyond this simple model by combining a measurement of the dark matter 
clustering, galaxy clustering, and their cross-correlation by defining a biasing b 
and cross-correlation r such that: 

where N ap is the galaxy number count density contrast smoothed with a com- 
pensated filter. Therefore, A ap is similar to M ap , except that it applies to the 
number count instead to the shear. As we discussed before, the compensated 
filter is a passband filter, quite narrow in Fourier space. If one chooses the num- 
ber count fluctuations -/V ap to be a foreground galaxy population with a narrow 
redshift distribution, then the biasing and cross-correlation b and r emerging 
from Ea(|55jl will be relatively localized in redshift AND wavelength. The com- 
binaison of well localized redshift and wavelength corresponds to a roughly fixed 
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physical distance. Therefore we can say that, even with the simple scheme of 
galaxy biasing given by Ea(|54|). an estimate of b and r from Ea (|55|) is fairly 
local in physical scale, for the foreground galaxy population under considera- 
tion (Schneider 1998, Van Waerbeke 1998). This result has been proved to be 
robust against a wide range of cosmological parameters and power spectra (Van 
Waerbeke 1998). 

This idea has been applied for the first time in the RCS data (Hoekstra et 
al. 2001). Unfortunately, this survey is not deep enough to provide an accu- 
rate measure of the dark matter clustering that could allow to separate b and r. 
Instead, the authors measured the ratio b/r = 1.05!to'io f° r the favored ACDM 
model (Slo = 0.3 and £l\ = 0.7). On the other hand, a combination of deep and 
shallow survey could help to measure the bias and the cross-correlation indepen- 
dently. This was done by combining the RCS and VIRMOS-DESCART surveys 
(Hoekstra et al., 2002). RCS is a wide shallow survey with a mean source red- 
shift of ~ 0.6, and VIRMOS-DESCART is a deep survey with a mean source 
redshift ~ 0.9. By selecting the foreground population with a median redshift 
~ 0.35 on the RCS survey, the number counts ), and the cross-correlation 
(M a pAap) were measured. The aperture mass (M^ p ) is measured on the deep 
survey. Figure E31 shows the measured b and r as a function of scale (angular 
scales are also converted to physical scale for a given cosmological model, with 
the lenses at z = 0.35). Although a proper interpretation of the measurement 
requires a better knowledge of the redshift distribution and cosmological pa- 
rameters, it is a direct indication of the stochasticity (r < 1) of the biasing at 
small scale, and that the biasing varies with scale as we approach the galactic 
scales, below 1'. The foreground galaxies were selected in R, and it was found 
that b = 0.7lig;^ on a scale 1 — 2 Ii^q Mpc, and r reaches a miminum value of 
r = 0.57ig;Q7, at 1 Mpc. We should note that b tends toward 1 at larger 
scale. 



6. Dark matter power spectrum inversion 

The central interest in cosmic shear observation is dark matter. This is probably 
even more important than measuring the cosmological parameters, for which we 
have some hope to measure them very accurately in the future (although there 
is the issue of degeneracies where lensing can help). One important question is 
then: what can we say about the dark matter distribution, provided we know 
all the cosmological parameters? This is nothing else but to try to map the 
dark matter in the same way we map the galaxies or the cosmic microwave 
background, or at least to measure its power spectrum in three dimensions, for 
all possible scales, independently of any evolution model. This is in principle 
possible from a direct inversion of Ea H25|) . but there are two issues here. One 
is that virtually, all physical wavelengths k are projected out to give a single 
angular wavelength s, and with a naive deprojection, one needs some cut-off 
somewhere in /c-space to perform the invertion. The other issue is that the 
3D power spectrum evolves non-linearly with redshift in the non-linear scales, 
therefore how could we be independent of any modeling when inverting the 2D 
convergence power? The first 2D convergence power spectrum estimate was 
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performed in Pen et al. (2002) on the VIRMOS-DESCART data, and in Brown 
et al. (2003) on the COMBO-17 data, but the spectrum inversion was not done. 

Pen et al. (2003) investigated the inversion using a singular decomposition 
technique, an extension of the minimum variance estimator deprojection devel- 
opped in Seljak (1998). The non-linear evolution of the 3-D power spectrum 
was assumed to evolve linearly with redshift even in the non-linear regime. This 
hypothesis is, surprisingly, a viable assumption within the scale range of interest, 
and produces errors still smaller than the statistical errors. The result is shown 
on Figure[2Ufor the dark matter (top) and the galaxies (bottom). It shows a very 
nice agreement with the cosmic microwave background C;'s (WMAP points ex- 
trapolated at z = 0, see Spergel et al. 2003), and with clustering measurements 
from other galaxy surveys. 

A dark matter-galaxy cross-correlation was also deprojected, allowing Pen 
et al. (2003) to estimate the 3D biasing b and matter-light correlation r. They 
found b = 1.33 ± 0.19 and r = 0.68 ± 0.24 for the /-selected galaxies. The 
bias value is slightly different than the one measured from the aperture mass 
on the RCS survey (section 5), but we should keep in mind that the galaxy 
populations are different (R compared to / selected galaxies for the RCS and 
VIRMOS surveys respectively). The physical scales probed in VIRMOS are also 
larger because it is a deeper survey than in RCS. 



7. Gravitational Lensing and Cosmic Microwave Background 

The use of lensing with other experiments improves the accuracy of cosmologi- 
cal parameter measurements and eventually breaks some intrinsic degeneracies 
attached to each. The potential interest of combining lensing by large scale 
structures and cosmic microwave background experiments has been studied in 
Hu & Tegmark (1999). The joint study of the weak lensing RCS survey and 
the WMAP data performed in Contaldi et al. (2003) is shown on Figure [2*51 and 
illustrates the gain of this combination: it provides a direct evidence of the low 
value of the matter density f2o 5 which indicates a high non-zero value for the 
cosmological constant, independently of the supernovae result. 



8. Approximations and Limitations 

Born approximation and lens-lens coupling The lensing theory developed in 
Section 1 assumes the lens can be projected onto a single plane, and therefore 
that the ray-tracing through a thick lens is equivalent to a thin lens appropriately 
weighted. As it has been quantified by Bernardeau et al. (1997), Schneider et 
al. (1998) or Van Waerbeke et al. (2002), it turns out to be a very good 
approximation. If we call 9 the direction of the unperturbed ray trajectory, a 
ray-light passing through a first lens will be slightly deflected by an angle SO, 
and will impact the second lens at a position angle 9 + 59 instead of 9 if the 
light ray were unperturbed. From a perturbative point of view, it means that 
expression Eq ()12j) has a correction term because the position angle to compute 
the lens strength is no longer x = fx(w)9, but 



Aussois's Cosmic Shear Lecture 



dark matter 



o.i 



0.01 



0.001 



0.01 





VIRMOS DESCART 

CMB 

RCS 



0.1 

k (h Mpc ') 



1000 



100 



10 



CO ^ 

< 



0.1 



0.01 



0.001 



0.01 



galaxy 



□ VIRMOS 
o 2MASS 
« APM 
• SDSS 



m ^ 

' I 

IE* |I | 







0.1 1 
k (h Mpc" 1 ) 



10 



Figure 24. Top: Dark matter 3D power spectrum, deprojected from 
the 2D convergence power spectrum measured in the VIRMOS survey, 
using SVD (Pen et al. 2003). The power is rescaled to z = 0, points 
are compared to the CMB (WMAP) and RCS lens survey. Bottom: 
Galaxy 3D power spectrum deprojected using the same method. For 
comparison, points from 2MASS, APM and SDSS are also shown. 
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Figure 25. The two dimensional, marginalized likelihoods for the 
(Q m ,as) plane. The overlaid, filled contours show the 68% and 95% 
integration levels for the distributions. Bottom - RCS only, Middle - 
CMB only, Top - CMB+RCS. Courtesy Contaldi et al. 2003. 
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Given that the correction to the light trajectory is a second order effect in the 
perturbation, it is expected to become important in any high order statistics of 
the lensing fields. Mathematically, indeed, they have the same order than the 
second order dynamical correction (which is proportional to the second order 
gravitational potential <I>( 2 )). It turns out that the light trajectory correction is 
much smaller than the dynamical second order correction. The reason is that 
Ea(|57jl involves a second lensing efficiency factor (the ratio of angular diameter 
distances /it's), which is not present in the second order dynamical correction. 

Figure shows several comparisons of the non-linear prediction for the sec- 
ond and third order statistics with a measurement of the same statistics done in 
ray-tracing simulations. It clearly demonstrates that non-linear calculations give 
quite accurate results, and that approximations related to the ray trajectories 
are valid to better than 2%. 
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Figure 26. Left panel: variance and skewness of the top-hat 
smoothed convergence field a OCDM (Oo = 0.3, as = 0.85, z s = 1.03, 
r = 0.21) and SCDM model (Q = 1., a 8 = 0.6, z s = 1.03, V = 0.5). 
The solid lines show the non-linear predictions, and the dashed lines 
the leading order of the perturbation theory calculations. The vertical 
dotted lines in the left panels denote the reliable scale limit fixed by 
the resolution of the ray-tracing simulation. The "large" error bars 
correspond to a survey of 25 square degrees, and the small error bars 
to 1000 square degrees. Right panel: Skewness correction due to 
the Born approximation and lens-lens coupling for SCDM (top) and 
OCDM1 (bottom) measured in ray-tracing simulation. 
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Non-linear lensing effects To first approximation, we consider the galaxy ellip- 
ticity an unbiased estimate of the shear. However, Ea (|43|) tells us that lensing 
is really non-linear. This approximation has been estimated in Barber (2002): 
it is negligible for sources at redshift less than z ~ 1 and for scales larger than 
5', while at smaller scale, a few percents effects could be detected. Fortunately, 
the use of the full non-linear lensing equation does not present any theoretical 
or technical difficulties, so small scale non linear lensing effect can be easily 
handled. It is just usually ignored in most of the theoretical and numerical 
works. 

Non-linear modeling On the other hand, Figure I2U also demonstrates that 
the accuracy of non-linear predictions on the 2-points statistics is never better 
better than 10%, while it is never better than ~ 20 — 30% for the skewness. 
This theoretical limit is a severe issue (Van Waerbeke et al. 2002) since one 
cannot expect to do precision cosmology if the accuracy of the model we use to 
extract the cosmological parameters is worse than the precision we want to reach 
on the cosmological parameters (that is a few percents). Smith et al. (2003) 
proposed an improved version of non-linear modeling, which is unfortunately 
still insufficient. In particular, to increase the precision, we still do not know 
whether the baryons have to be taken into account in the modeling or not. The 
goal here is a modeling accurate to 1 — 3%, if one wants to reach the same 
accuracy on the cosmological parameters. 

Intrinsic alignment Gravitational lensing is not the only natural process which 
produces alignment of galaxies over large distances. Intrinsic alignment might 
occur from tidal fields, and produce galaxy shape correlations over cosmological 
distances, and contaminate cosmological signal (Croft & Metzler 2001, Catelan 
k Porciani 2001, Heavens et al. 2000, Catelan et al. 2001, Hatton & Ninin 
2001) which should in principle split, in a predictable way, into E and B modes. 
There is unfortunately only partial agreement between the different predictions. 
Moreover, most of the predictions stand for dark matter halos, while we are in 
fact observing galaxies, which should experience some alignment mixing. This 
has not been simulated so far. Concerning the dark matter halos alignment, 
despite the disagreement among the predictions, it is generally not believed to 
be higher than a 10% contamination for a lensing survey with a mean source 
redshift at z s = 1. An exception is Jing (2002), who suggested that intrinsic 
alignment could dominate the cosmic shear even in deep surveys. This possibility 
is already ruled out by observations: this would indeed imply a very low erg ~ 0.1 
if the observed signal were dominated by intrinsic alignment, and we should also 
observe an increase of the effect as we go from deep to shallow survey, which 
is not the case (see Figure IT7j) . In any case, intrinsic alignment contamination 
might be an issue for studies using a single source redshift in their analysis. In the 
future, this will not be the case since photometric redshifts will be available. In 
that case, the effect can be suppressed by measuring the cosmic shear correlation 
between distant redshift sources, instead of measuring the fully projected signal. 
Consequently, intrinsic alignment should not be considered as a critical issue 
(Heymans & Heavens 2003, King & Schneider 2003). Pen et al. (2000) and 
Brown et al. (2002) reported the first two evidences for intrinsic alignment in 
the nearby Universe, which are not too inconsistent with the predictions. 
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Source clustering Source clustering arises because a subset of sources overlap 
with a subset of the lenses which are probed. There is therefore a natural bias to 
measure the signal preferentially in high density regions, across the overlap area. 
This effect gives rise to correction terms in high order statistics (Bernardeau 
1998). It is easy to understand the problem if we model the source redshift 
distribution including a clustering term: 

p w (e,w) = p w (w)(l + 5 ga i(f K (w)e,w) + ...), (58) 

which replace the source redshift distribution p w (w) in Ea(|17|). It is then easy 
to see that a density coupling occurs in Eq ()16l) . The source clustering effect 
was extensively studied by Hamana et al. (2002). They confirmed that it is not 
an issue for the 2-points statistics, but could be as high as 10 — 20% for the 
skewness of the convergence, for a narrow redshift distribution. In case of the 
broad redshift distribution, the effect is diluted by the bulk of non-overlapping 
areas. For future surveys, an accurate measure of the high order statistics will 
require a precise estimation of this effect, which is not a problem by itself, 
but it must be done. Schneider et al. (2003) predicted that source clustering 
could produce B mode in the shear maps, which has not being tested against 
reay-tracing simulation yet. However, the predicted amplitude is 2 orders of 
magnitude and below for aperture mass scales larger than 1'. 

PSF correction With the non-linear modeling of the power at small scale, this 
is certainly the most serious issue concerning the cosmological interpretation of 
the cosmic shear signal. Again, if we want to reach a few percents accuracy 
on cosmological parameters measurements, we need a PSF correction with that 
accuracy. So far we are able to reach 10% precision with the KSB method for 
a typical signal measured on sources at z = 1 (Erben et al. 2001, Bacon et al. 
2001, Van Waerbeke et al. 2002). This is reasonably good, but we still need 
to gain one order of magnitude (in addition to the order of magnitude we need 
to gain for the non-linear modeling as well). The 10% uncertainty is an upper 
limit, which comes from the large B mode found in all surveys, at different 
scales, for probably different reasons (for instance, RCS have B mode at small 
scale only they may have measured intrinsic alignment). This upper limit is 
reduced if one uses the scales with very small or no B mode, but then some 
cosmological information is lost. So far, our understanding of the PSF modeling 
is insufficient in particular concerning the PSF variation (and stability) accros 
the CCD's and the contribution of high frequency modes. Space data are often 
viewed as potentially systematics-free. This is unfortunately not true, since all 
space data which have been processed for cosmic shear, required a significant 
PSF correction. However, the main difference between space and ground based 
data is that, in space, the PSF is certainly more stable between exposures. But 
one should not forget that in space, the PSF is 100% instrumental (it is the Airy 
spot, which is larger than the Airy spot on the ground because space telescopes 
are smaller), and not atmospheric at all (which it is with ground based data 
with larger telescopes). Dealing with a non circular Airy spot to correct for 
the galaxy shapes was not trivial for the Hubble Space Telescope for instance, 
mainly because of the severe undersampling of the PSF (Hoekstra et al. 1998). 
There is no intensive simulation of shear measurement under various realistic 
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Figure 27. Top hat variance of shear as function of angular scale 
from 6 cosmic shear surveys. The open black stars are the predictions 
for the CFHT-LS which will start by 2003 with Megacam at CFHT. 
This is the expected signal from the "Wide Survey" which will cover 
170 deg 2 up to Iab = 24.5. For most points the errors are smaller than 
the stars. 



space image conditions, only qualitative estimations have been done (Refregier 
et al. 2003, Massey et al. 2003), which seems promising. 

Finally, one should emphasize that the most difficult part of the PSF cor- 
rection is not the anisotropic correction, which is done quite accurately, but the 
isotropic correction (Erben et al. 2001, Hirata & Seljak 2003). The ultimate 
limit of PSF correction in space and on the ground is still an open question. 



9. Prospects 

In the WMAP (Spergel et al. 2003) context, one can wonder whether future cos- 
mic shear surveys can still provide useful cosmological informations that would 
not be available otherwise from CMB and SNIa experiments. The answer is 
clearly yes because cosmic shear is the only way to directly probe dark matter 
on scales that cannot be directly probed by other techniques. It can explore 
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Figure 28. Theoretical expectations on cosmological models beyond 
the standard model from the wide CFHT Legacy Survey that will cover 
170 deg 2 . The dots with error bars are the expected measurements of 
cosmic shear CFHTLS data. The lines shows various models discussed 
by Benabed & Bernardeau (2001). 



the properties of the dark matter together with luminous matter on quasi-linear 
and non-linear scales where the complexity of physical processes make theoret- 
ical and numerical predictions among the most challenging tasks for the next 
decade. Further, because cosmic shear is sensitive to the growth rate of per- 
turbations integrated along the line-of-sight, the additional redshift information 
provides a tool to study the structure formation mechanism and the clustering 
history with look-back time. It is the purpose of tomography to study the 3D 
matter distribution by combining the lensing effect with the redshift information 
of the sources (Hu 1999, Heavens 2003). This clearly belongs to the prospective 
part, and was not discussed in the review. 

The future key scientific goals are therefore the reconstruction of the 3- 
dimension dark matter power spectrum as function of redshift, the analysis of 
the properties of the relation between light and mass, and the study of the dark 
energy equation of state. As shown in this review, the scientific studies of the 
cosmic shear surveys done so far already show that there are neither conceptual 
nor technical barriers that hamper these goals to be achieved quickly. Tegmark 
&: Zaldariaga (2002) with the RCS cosmic shear survey and Pen et al. (2003) 
with the VIRMOS-DESCART cosmic shear surveys have demonstrated that the 
3-D power spectrum of the dark matter can be reconstructed. Remarkably, 
their results extend monotonically toward small scales the dark matter power 
spectrum derived from CMB experiments. Likewise, Hoekstra et al. (2001) 
Hoekstra et al. (2002), and Pen et al. (2003), have shown that the properties of 
the biasing and the dark matter-galaxy cross correlation can already be analyzed 
with present-day surveys. Finally, Bernardeau et al. (2003) as well as Pen et 
al. (2003) have shown that high order statistics are already measurable from 
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ground based data covering 10 deg 2 , thus providing independent informations 
on cosmological models, with eventually some important degeneracies broken. 

Although the B-mode contamination is still an important technical issue 
that may slow down the cosmic shear developments, in principle the next large, 
deep and multi-color surveys will be in position to address questions relevant 
for cosmology and fundamental physics with a high degree of precision. Such 
surveys, covering hundreds of degrees, with multi-bands data are about to start. 
Among those, the wide CFHT Legacy Survey 5 will cover 170 deg 2 , spread over 
three uncorrelated fields, in 5 optical bands, and a fraction will be followed 
up later in J and K bands with the wide field infrared camera WIRCAM at 
CFHT. Figures H3 and H3 show some predictions of CFHTLS. On Figure 1771 
we simulated the expected signal to noise of the shear variance as function of 
angular scale for a ACDM cosmology. The error bars are much smaller than 
the VIRMOS-DESCART survey which has the same depth as CFHTLS. On 
Figure 1281 we compare the expectations of the CFHT Legacy Survey angular 
power spectrum with the predictions of several theoretical quintessence fields 
models. It shows that 200 deg 2 deep survey with multi-color informations to get 
redshift of sources, one can already interpret cosmological data beyond standard 
interpretations. The CFHTLS will be of considerable interest because one of the 
fields is also a target for the VMOS/VDDS spectroscopic survey (Le Fevre et al 
2003), the XMM-LSS survey (Pierre et al 2001) and also the COSMOS Treasury 
Survey that will be done by the HST/ACS instrument. Hence, in addition to a 
complete description of the redshift distribution of the CFHTLS galaxies, as well 
as of the X-ray clusters and active galaxies, high accuracy shape measurement 
of galaxies will be feasible. This HST/ACS data set attached to a subsample of 
CFHTLS data will permit to check the reliability of ground based PSF corrected 
shear catalogs but also to extend the shear analysis on very small scales, down 
to the galactic dark halos scales. Join together with CMB, we then expect to 
get by 2005 a complete view of the dark matter power spectrum and the biasing 
from Gigaparsecs down to kiloparsecs scales, as well as a detailed description 
of individual dark halo properties and of the redshift distribution of lenses and 
sources (Cooray & Sheth, 2002). One should mention the extensive search for 
dark matter halos in the Suprime-Cam fields (Miyazaki et al. 2003): this is more 
a cluster physics-oriented, but it illustrates very well the future developments 
with cosmic shear data for understanding the halo distribution and formation. 

CFHTLS is one of the new generation surveys, with similar studies begining 
at SUBARU, soon at ESO, with the VST, later with VISTA, also the NOAO 
deep survey 6 , Dark Matter Telescope 7 , and the PAN-STARRS 8 . Beyond 2005, 
space based dark energy/matter probes like SNAP appear as a kind of final 
achievement. In principle, SNAP can provide deep images, accurate photometric 
redshift, a large field of view and outstanding image quality one expect for 
cosmic shear. A dark matter space telescope, entirely dedicated to cosmic shear 



5 http://wow.c£ht.Hawaii.edu/Science/CFHLS/ 
6 http: //www. noao.edu/noao/noaodeep/ 
7 http://www. dmtelescope.org/dark_home. html 
8 http: //www. ifa.hawaii.edu/pan-starrs/ 
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observations, might also be an interesting option: it could be a 'small' telescope 
(therefore fairly cheap) that could observe the shear over all the sky. One could 
then 'see' the dark matter everywhere! 

Cosmic shear data are optimized when they are used together with other 
surveys, like Boomerang, CBI, DASI, WMAP of Planck CMB experiments, SNIa 
surveys, or galaxy surveys (2dF, SDSS). The first tentative recently done by 
Contaldi, Hoekstra & Lewis (2003) shows that tight constraints can really be 
expected in the future. Likewise, by using cosmic magnification instead of cosmic 
shear on the 100, 000 SDSS quasars, Menard & Bartelmann (2002) have shown 
the cross-correlations between the foreground galaxy distribution and the quasar 
sample is also useful to explore the properties of the biasing. In principle magni- 
fication bias in the SDSS quasar sample can provide similar constrains as cosmic 
shear. Yet, this is a widely unexplored road. 

This review shows that the Aussois winter school has been organized at a key 
period, between the first and second generations of cosmic shear surveys. The 
first period, from 1999 to 2003, demonstrated cosmic shear can be detected and 
exploited for a lot of cosmological questions. We now enter the next generation 
surveys, which will end around 2010. These are large ground based surveys (like 
the CFHTLS-weak lensing) which will fully exploit the new windows opened by 
the first surveys. Then we will enter the third (last?) period with extensive 
space observations, probably after 2010, like SNAP, which will permit to do 
precision cosmology, and maybe to close the subject. 
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